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TiiGre is airailabl© a program v^hieh will so1%''g linear? pro^- 
granming px'oblen^ by the simplex method, 'Biis program has beem used 
with good sue ess 3 on problem 216 « 

DeseriDtion of ProblBia 

This program, is used to find the imicisrom of a linear func- 
tion in several Yariables, 
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subject to the restrictions 
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and x^ - Oj j - 1, — ^ m, 



This differs frora the more general forimilation of the problem in that 
no 0q.-aalities are al3.07#sa imder (S)« It is also required that the 
initial irariable values be zeros « 

Us© of the Program 



The program starts with the matr^^ix of coeffieieatSj (a.^. .J j 
stored by columae oa the auxiliarij drum beginning with as. arbi« ~*'^ 
trar^^- drura afldrass^ vghieh is the preset parajieter ph 1» Tub "value 
of the fiiactioa^ initially zero, is stored in CM at flad 21 ^ which 
isiist be assigned « ITollowing it,, starting at SZlg, are the zmmbsrs 
-e.j i^e^g ths negatives of the eosffieisnts in the linear foxTti (!)„ 
TM4. 3?estriction eoEstaats b. are also in (M^ startii3.g ?d.th flad mlj 
whioh iwdBt be assigned* la addition ^ stoi^age for two coliHime of 
the siatrix must ba proYided at flacls pi ami mSj both of which must 
be assigned e All of these xivMhers are la (g^^S) foi'm* 

Besides these ^ the pTOgi-am recjulres a block of WI iategero 
to k©sp track of ths order of the variables and rest riet ions « There 
Yiill bs ii^s of these single register umabes^e^ referred to as tags^ . 
and it is' arbitrary? 'i^hat mmbers are ussdj although it is Tsry G'Ob- 
Yenient to us© :eela.t±Ye addresses of the Tariablss is stibsequent prO'- 
grams* The tags for th© restrictions start at flad dig which must be 






assigEfsdj, and lieBcs',. tlie tags for the variables start at dl-^-ji., 

Finall^^ the para^riDterei pi-Ll •-•" n ~ nmnber of variables. 
p.ml ~ m =- niMber of b., *s and pfel = chraiii. address of niati^jji' must be 
assigned. A2.teraat8ly . these laaj^ .be isieerted in three reglBters of 
the program.,, m at SOr, phi at Sir., aad n at SSr^ 

TixG entry to tha prog3?ain is a ta instruction at os% The 
exit is an sp to the original control miless an infinito s^xj-mim ±b 
indicated J in x'^iiich ease a jujip to is made from 4t8» 

'FoT-m of Results 

The niaxisBl Taluo of s will bs at si. 

The values of the Tariables, x.. , and the amoiirits by whicli 
the restrictions (2) are inaqualitiaSj i'ie^j, the so^'CsJ.lad jils^^k 
jJBriab^leSa p^. - b.. - Sl^'-.- .• ^~-f ^^■■•^- d.et®Tm±ned. by the qtiaBtitios left 
at M^'sjii by^the order of^thd tags* Recall that the tags corrospond 
in one-to-one fashion with the -variables x, and slack variableB p.* 
To find the desired values the foILloyjiag cfo^aspondenee is ussdi ™ 

(milder at ml -'■■ £1)4—=^--^ (variable or slack -rariable 

eorresDondlBE to tag at fll-^i) 
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■ (ssro) itj»«_.«.„.„,„»»^»„^Hg' ("^as^'iable or eUack -UBriabla 

correspoMing to tag at dl-^iH-j) 

if aMrssBos (i^elatiTo or othsrwiss) are usaci as tags, a sorting 
techniqiia which picks up the a sioiibers at ml aiKl places tliem at tho 
addresses ^ich are the correspoiiding tags and thou fills in the 
reniainiiig aclclressos with s^eros is coiiueaiosit^ 

Other results which are soi!istiio.as ussfi are the iimribsrs 
vfhioh Charaes denotes by s. - c.* Siese will bo found starting: at 
Szl in an order corresponding to tha last in tags a, 

'Ehe tiJiBX iBStrix is, the Dsrtioily ia¥©3:?te5. os^iginal obSj. 
and the order of rows and cojAiiiins also eorresponds to the final orcles? 

Oi.' til© 'C8.gSo' 

For coaTeniencOs the niMbor of passes is counted in both 
FF6 and regis tex" tl3^' In additioiij tha ¥-a3.u© of s at aiij moiTieat of 
the ffiaximisiag procase la stored as a {243.6} nmsber in FF4 and FFS^ 
Bins thes© results igay fee obtained as the process coatimiosaM can 
be usoci to dstesmiae a stopping point previous to the attaiaing of 
the actual masiKms* la ordsr that the program eouM ho s topped at 
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operating Time 



The operating time is deteriiined by tira factors, one of 
lAtiiicli can be iiisasured preeisolsr, l^e^^ the time per passj and the 
other of which can orO^y be estiimted^ ieG«9 th© nmnb&r of passes* 
Since the time consmiied in a pass is mostly due to pe3?foriiiing an 
algorithm on each element of aa iixm rmtrlx, the time should he 
proportional to mIo Upon actually plotting pass times vs« mn a 
suitable equation was found; 

tine per pass « «000075ina •*• .03 minutes. 

'Bie classical estimate of the number of passes is 2m« In Involve 
different trials t^ith niatxuces ranging from 31 x 10 to 77 s: 36 the 
number of passes was foimd to vary from 1.2m to 2,4m<j with an 
apparent tendency to give the fevier number of passes vdth smaller 
matrices. In addition, t?mnty~five trials ¥%'9r8 made in v^hich the 
initial conditions were somevifhat closer to the maziimni, with the 
nuBtber of passes varying from lelEi to l.Sm. 

Notes 



If it is desired to eliMnate the restriction z. §, for 
one or inore variables the classical procedure has been *^ to make a 
variable substitution which doubles the number of variables* In 
problem. S16 it vms found convenient to group the variables accord* 
ing to order of njagnitude and then make on3-y one substitution for 
each group. For instance, if say^ x,. » .*• x^ fiere to be noa« 
restricted and {x^ , Xo),? (x„ * x.^ Xgt were tne grouping according 
to order of inagni-sude*^ then tv;o ne?/ x?ariables y- and y^ Yjould be 
introduced, and the variables used in the program would be 

^1 "^"' ^1» ^S "^ ^1' ^3 "^ ^2* ^4 '^ "^2* ^5 '^^g* ^^1* ^^2 

©lese seven variables ¥irould be restricted, i.e.y x^ •«- f^t 0^ 
but the restrictions x^ ^. -y., , 8tc*g will alioii 



x^9 *.«, Xg to become as negative Us necessary < 

Degeneracy has heen said to be the greatest fault of the 
siniplex method 9 to the possibility 'of ^cycling^*.* Howeverj this 
possibility is quite r-are, so no precaution against it was imA® 
in the program* It was fcund^ hoimverj that the closely comiae- 
ted question of linear independence of variab3.es could give 
trouble 5 a false infinite mazimua caused by round-off error*, It 
was found that in the case vihere several variables were linearly 
dependent 5 as in the escample above ^ that due to round^off error 
they did not appear to be so 5. and thus 'could occur in the basis 
simultaneously a To see how this could cause a false infinite maxi^ 
rnxm consider the esaii5)le above*, Suppose that z- -^ y»j x^ -^ y^ ^ 
and jj vmre all in the basis, althoug:h they are linearly dependent, 
KoVif tne coefficients of these variables in the linear form s are 
c^ , Cgg. and -e^'-Cg respectively^ If^ however, -e^-Cp had acquired 
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excess, e, of round-off error, xm could increase y without 
md causing z to increas© also; i*©*,, if y^ == M, a large 
nujribei', th©n the computed s v/ould look like 

z' « c- (^-^-M) •»• CgCXg-^M) -^ ("Cg-c^-J-e)!!*! -^ other t©ims 

« c- Xj -5- Cg Xg •«- Me + other terms ,, 

in Spite of the fact that x- and Xg are such that restrictions (S) 
are satisfied • Since this raise uiaziiiiuia actually occurred ^ a pre- 
caution should be taken against it. It is sufficient to decrease 
the coefficient of y. in the form (1) by an aiaount larger than 
any e3cpected round-off error but not enough to seriously effect 
the value of z* 
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